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ABSTRACT 

We present high resolution imaging of the puzzhng radio and optical nebula 
G70.7+1.2 with the Keck Observatory's laser guide star adaptive optics (LGS- 
AO) system and the Chandra X-ray Observatory. The archival X-ray observations 
show a hard (F ^ 1.8), low luminosity [Lx ~ 4 x 10^^ erg s~^) point source at the 
center of the nebula. Follow-up LGS-AO near-infrared imaging of the Chandra 
error circle reveals a relatively bright {K' ?5i 14 magnitude) counterpart. Both 
its color and brightness are consistent with a heavily obscured B-star or possibly 
a late-G/early-K giant. The most plausible explanation is that this newly dis- 
covered X-ray source is a non-accreting B-star/pulsar binary powering the radio 
and optical nebula. If so, the luminous Be-star discussed in the literature seem- 
ingly embedded in the nebula is not the dominant force responsible for shaping 
G70.7+1.2. Thus, we suggest that G70.7-I-1.2 is the result of two unrelated ob- 
jects (a B-star X-ray binary and a Be star) interacting with a dense molecular 
cloud. With this explanation we believe we have solved the mystery of the origin 
of G70.7+1.2. 

Subject headings: X-rays: binaries — ISM: individual (G70.7+1.2) 



Introduction 



G70.7+1.2 is a cor npact 20") radi o and optical nebula in the Galactic Plane whose 



origin is controversial ( IReich et al.l Il985l : iGreenI Il986l : Ide Muizon et al.l Il988l : iBally et al. 



19891 ). The shell-like radio morphology is accompani ed by broad, blue-shifted [O I] and [S 



II] em ission lines indicative of an interstellar shock (Ide Muizon et al.lll988l : iKulkarni et al. 



19921 ). Millimeter CO emission traces t his radio and opti c al emission, suggest i ng the shock 



is int eracting with molecular material (IBally et al.l Il989l : iPhillips et al.l Il993l : lOnello et al. 



19951 ). In addition, a bright near-infrared (NIR) star appears to be embedded in G70.7-I-1.2, 
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and it is surrounded by a strong Ha reflection of spectral type Be (IBecker fc Fesenlll988 
Kulkarni et al1ll992h . 



Ironically, it is the plethora of clues t hat make G70.7+ 1.2 a perplexing object, despite 
its discovery more than two decades ago (IReich et al.l 119851 ) . G70.7+1.2 is unique in that 
nearly every Galactic prototype has been proposed to explain it: young s upernova rernnant 



nova shell, stellar wind bubble, H I I region and Herbig-Ha ro-like outflow (iReich et al.lll985 



de Muizon et al.lll988l : lGreenlll986l : iBecker fc Feserull988l ). However, none of these can ex- 
plain the low expansion velocities and the non-thermal radio emission. The only currently 
proposed consistent theory to explain these properties is one in which the bright NIR star is 
paired with an unseen r ieutron star to form a Be-radio pulsar binary moving supersonically 
through the dense gas (IKulkarni et al.lll992l ). In this model, significant mass loss from the 
luminous Be-star inflates a bubble which is filled by a mixture of the stellar wind with ener- 
getic particles and magnetic field from the pulsar. This combination creates the non-thermal 
radio emission coincident with the optical bow shock of the medium surrounding the system. 
This model makes the prediction that an X-ray source or pulsar should be seen coincident 
with the embedded Be-star. 

Here, we report on archival X-ray and new Keck LGS-AO observations that reveal an 
X-ray source with a NIR counterpart in the center of G70.7+1.2 which is unassociated with 
the Be-star. We suggest that the object known as G70.7+1.2 is the result of the interaction 
of a luminous Be-star and an X-ray emitting B-star/pulsar binary with dense molecular 
material. If correct the resulting study of this object will be an interesting laboratory for the 
study of plasma processes. The observations and results are presented in fj2l In ^we discuss 
the implications of this source as it relates to resolving the mystery surrounding G70.7-I-1.2. 



2. Observations and Analysis 
2.1. X-ray 

G70.7+1.2 was observed 2003 October 11.33 UT with the ACIS-S detector on Chandra 
in the standard, timed exposure mode. The archival data were analyzed with CIAO version 
3.^. We reprocessed the level 1 events from the Chandra X-ray Center (CXC) in order to 
make use of the latest calibration and removed pixel randomization. The level 2 event file 
was created by filtering grades 0,2,3,4,6 and good-time intervals. The total exposure time 
after filtering periods higher than 3-0" above the mean background level was 37.6 ksec. 



http: / / www.cxc.harvard.edu/ciao/ 
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Diffuse emission and a point source (hereafter CX O J200423.4+333907; Fig ure [T]) are 
detected at the position of G70.7+1.2 (as first noted by lArzoumanian et al.ll2004l ). We com- 
pared the positions of 27 X-ray sources on the S3 ch ip with counterpa rts in the 2MASS 
point source catalog to correct the native astrometry (jCutri et al.ll2003l ). This comparison 
showed evidence for a small systematic shift, Aq;2mass-cxo = -0'.'13 ± (X'll, A52mass-cxo = 
-O'.'OS ± O'.'ll. The best-fit position of CXO J200423.4+333907 including this offset is 
a(J2000)=20'^04°^23.430'^ and 5(J2000)=33°39'06'/73 with measurement uncertainty of 0'.'03 
and O'.'OS in each coordinate, respectively. Combining the measurement and transformation 
errors in quadrature gives an uncertainty of O'.'IS (l-cr) for the X-ray position of CXO J200423. 4+333907 
relative to 2MASS. This po sition lies 3'.'6 from the nearby luminous NIR st ar, which implies 
the two are not associated (IKulkarni et al.lll992l : lArzoumanian et al.l 120041 ). 



We applied the adaptive smoothing algorithm CSMOOTH to highlight the diffuse emis- 
sion after subtraction of CXO J200423. 4+333907 and produced a fiux-calibrated image by 
applying an exposure map as outlined in the CIAO threads. The contours of this emission 
are overlaid on an NIR image (see §2.2p of the nebula in Figure [1] Evidently, most of the 
diffuse X-ray emission is not coincident with the diffuse NIR emission. 

We extracted photons within a I'.'S circle (corresponding to 90% of the expected counts 
at 1.4 keV) around CXO J200423. 4+333907 to perform spectral and variability analyses. 
The source contains only •i^'^^l counts. Upon examination of regions both inside the diffuse 
emission and in a source-free area, we expect only two of these to be background photons. 
A Kolmogorov-Smirnov test shows that the arrival times of the source photons differed from 
a constant rate at only the ~ l-cx level, thus the source cannot be considered variable. 

We begin our spectral analysis by noting that all the photons from CXO J200423. 4+333907 
fall in the range 1.0-4.5 keV, suggesting a hard spectrum. After calculating the response ma- 
trix and effective area of this portion of the CCD, we fit an absorbed po wer-law mo del to 
the spectrum using Cash statistics (due to the limited number of counts; ICashlll979l ). The 
best-fit parameters in Tabled] show a relatively hard photon index, F f» 1.8, an d low lumi- 
nosit y, Lx ~ 4 X 10'^^ erg s'^ (2.0-10.0 keV) for an assumed distance of 4.5 kpc (IBally et al. 



19891). These values are consistent w i th kn own accreting neutron stars in quiescence (e.g. 
Rutledge etaPboOll : ICampana et~aDl2005h . In addition, we fit two absorbed power-law 
models with fixed parameters (see Tabled]). The first has the photon index set to a typical 
value for quiescent neutron stars, F = 2, while the second has the column density fixed to 
the best fit value of the diffuse emission, A'^^^ = 1.0 x lO^^cm"^ (see below). 

The probability of finding a source as bright or brighter than CXO J200423. 4+333907 
within the extent of G70.7+1.2 can be determined from the local source density. A WAVDET 
analysis of the active CCDs (ACIS-12,3 and ACIS-S1,2,3,4) finds that 12 sources are as bright 
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or brighter than CXO J200423. 4+333907. The inferred density is then ^112 sou rces /deg 



This density is consistent with observations taken as part of the ChaMPlane Survey (iGrindlay et al. 
20051). which predi cts ~ 100 sources/deg^ with fluxes as bright or brighter than CXO J200423. 4+333907 



(IHong et al.ll2005l ). Consequently, there is a 0.3% probability that such a source would be 
found within G70.7+1.2 by chance. 

The diffuse emission presented enough counts for basic spectroscopy with statistics. 
We extracted events from a region of dimension ^ 30"x30" surrounding the diffuse emission 
(excluding the point source) and a source-free background region immediately east of the 
nebula with the same shape. This yielded 690+26 source counts, of which ^ 320 are expected 
to be due to the background. The resulting source plus background photons were grouped 
such that each bin contained at least 25 counts. 

The background subtracted spectrum was analyzed using XSPECvll^. We fit two 
models modified by absorption to the spectrum: a power-law and a Raymond-Smith plasma 
(see Table [1]). The unphysically steep photon index of the power-law model and the lower 
xt value lead us to adopt the Raymond-Smith model for the remainder of our analysis. 
The derived value of Nh is reasonably c onsistent with that of CXO J200423. 4+333907 and 
estimated value of 1.25 x lO^^cm"^ from lDickey &: Lockmaru (jl990l ). Integrating this model 



over the 0.5-2.5 keV bandpass implies a luminosity of 8.1 x 10^^ erg s ^ at the distance of 
G70.7+1.2, although it is not clear that this emission is associated with the nebula (see 



2.2. Near-Infrared LGS-AO 



G70.7+1.2 was observed under ph otometric conditions on 2005 April 30 UT wit h Laser 
Guide Star Adaptive Optics (LGS-AO: IWizinowich et al.ll2006l : Ivan Dam et al.ll2006l ) on the 
Keck II telescope and the Near-Infrared Camera 2 (NIRC2). We imaged the field in the J, 
H and i^'-bands with the wide camera of NIRC2, which provides a f«40"x40" field of view 
and a ~ 0'.'04 pixel scale. The H and /T'-band data sets consisted of five frames in each 
band. Each frame was exposed for 5 sec with 10 additions performed on the chip at five 
dither positions separated by ~ 30". The J-band data consisted of two images at the center 
of the chip. 

Each frame was fiat-fielded, background subtracted, and repaired for bad pixels using 
custom PyRAF softwar^. We then performed a second round of sky subtraction using a 



^ See http : / /heasarc.nasa. gov /docs/ xanadu / xspec / 

^PyRAF is a product of Space Telescope Science Institute, which is operated by AURA for NASA. 
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median combination of similarly processed frames of a nearby field. We used these processed 
images of G70.7+1.2 for photometric analysis, but produced a separate set of images for 
astrometry due to optical distortion in the NIRC2 camera. The distortion in the second set 
was corrected using algorithms derived from the preshipment review document^ with the 
IDL procedure provided by the Keck Observatorjo. The correction does not conserve flux, 
and thus is not suitable for photometry. 

We registered a median combination of the distortion corrected if-band frames to the 
2MASS point source catalog using 8 stars that were not over-exposed. We find residuals of 
0'.'04 and (/.'OQ in right ascension and declination, respectively. Registering the J and K'- 
band frames to this image yielded negligible residuals. Combining these errors with those 
in the X-ray position of CXO J200423.4+333907 yields an uncertainty of (/.'IQ (1-a) of the 
X-ray image with respect to the NIR images. Figure [T] shows the registered if-band frame 
with the Chandra error circle (99% confidence). We clearly identify a single bright NIR 
counterpart in all filters within the X-ray error circle. The best fit position of this source 
is a(J2000)=20^^04'"23.446^ and 5(J2000)=33°39'06'.'62 with an uncertainty of 0'.'04 and 0'.'09 
(relative to 2MASS), respectively. The centroiding errors are negligible. This position lies 
0'.'23 from the Chandra position. 

We performed aperture photometry of the source in each band on each individual frame 
relative to 2MASS stars in the field with the IRAF package APPHOT. We assume that 
the color term used to transform from the 2MASS Kg filter to the K' is negligible for our 
purposes. The uncertainties were determined with the 2MASS photometric uncertainty, the 
standard deviation of the zero-point determinations from the same 2MASS star in multiple 
frames and the photometric error of the NIR source itself added in quadrature. We find 
magnitudes of J = 15.56 ± 0.09, H = 14.51 ± 0.11 and K' = 13.97 ± 0.06. 

The probability of finding a star with K' ^ 14 magnitude in our Chandra error circle by 
chance is very low. To quantify this we extracted all sources present in the 2MASS catalog 
within 20' of G70.7+1.2. We find that the differential number of sources per magnitude per 
square arcsecond is well described by a single power-law with index 0.35 over the magnitude 
range 3 < K^, < 15. W e can conservatively assume (based on Galactic star count models by 



Nakajima et al.ll2000l ) that this can be extrapolated to our 5-0" detection limit of mx' ~ 20.0 
magnitude. From this we calculate that there is a ~ 25% percent chance of finding a source 
brighter than our detection limit in a circular region with a 0'.'49 (99% confidence) radius. 
However, the probability of finding a source with Kg = 14.0 magnitude or brighter is <^ 0.1%. 



"^available at http://www2.keck.hawaii.edu/inst/nirc2/ 
^See http://www2.keck.hawaii.edu/optics/lgsao/software/ 
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Thus it is unlikely that our NIR counterpart is drawn from the background population, and 
we assume that it associated with CXO J200423. 4+333907. 

The key issue in determining the nat ure of this source is the assumed extinction. Based 



on the colors of the luminous NIR star, iBecker &: FesenI (119881 ) estimate Ay ~ 5.6. This 



agrees well with the value of Ay ~ 5.4 obtained by taking Nh as d etermined from the 



spect rum of the diffuse X-ray emission and translating it into extinction ( jPredehl fc Schmitt 



19951 ). If we plot the NIR counterpart on a color- magnitude diagram (see Figure [2]) using 
this reddening we find that the star is consistent with a late G/early K giant spectral type 
at a distance of ~ llkpc. Consequently, the star is under luminous if it is associated with 
G70.7-I-1.2 at a distance of 4.5 kpc. 

The spectral fitting of the point source spectrum itself, albeit with poor statistics, 
implies a higher extinction of Ay = SA'^f'^. This allows for the possibility that the NIR 
counterpart is a heavily obscured main sequence B-star with Ay ~ 10.0 at the distance of 
G70.7+1.2. We prefer this interpretation when we consider the probabilistic arguments and 
existing multifrequency observations of G 70. 7+ 1.2 (see 



3. Discussion and Conclusions 



We identify a low-luminosity, hard X-ray point source with a NIR counterpart at the 
center of G70.7+1.2 using high resolution imaging. Both the measured X-ray luminosity, 
Lx ~ 4 X 10^1 erg s-\ and the photon index, F ^ 1.8, of CXO J2004 2 3.4+333907 are 



consi stent with quiescent neutron star systems (e.g. iRutledge et al.ll200ll : ICampana et al. 



20051 ). The magnitude and J-K' color in combination with the X-ray column density sug- 
gests the NIR counterpart is either an evolved background star or a heavily extincted B- 
star. However, an isolated background late-G/early-K giant cannot explain the observed 
X-ray flux. These stars have de ep convectiy e zones that power coronal X-ray emission, 
but it is typically ^ 10^^ erg s-^ flGiidell 120041 ). This is an order of magnitude below the 
required Lx ~ 2.5 x 10^^ erg s~^calculated assuming the observed X-ray flux at a dis- 
tance of llkpc. In addi tion, spectral types lat er than B2 have have observed X-ray lu- 
minosities ^ 10^^ erg s~^( lBerghoefer et al.l 119971 ). This suggests that the NIR source and 
CXO J200423. 4+333907 constitute an X-ray binary, and probabilistic arguments suggest 
that this binary is associated with G70.7+1.2. 

A simple geometric model can explain the existing multifrequency data (Figure [3]). The 
velocity of the molecul ar gas as measure d by CO observations is 5 km s~^ with respect to the 
local standard of rest (IBally et al.lll989l ). The stellar Ha line profile from the bright Be-star 
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is redshifted with respect to the CO with a velocity of 20-60 km s~ , while Ha reflected by 
dust in the eastern region is also redshifted with respect to the CO, but is blueward of the 
stellar Ha by 10-50 km s~^. This implies that the bright NIR star is moving into the cloud, 
away from the Earth. However, the [O I] and Ha throughout the rest of nebula traces the 
non-thermal radio emission and is uniformly blue-shifted by 10-120 km s~^ with respect to 
CO, suggesting that the source responsible for the shock is moving into the cloud, toward 
the Earth. 

The cloud size, as inferred from CO, is 3/^4,5 pc on the sky, where D4.5 is the distance 
to G70.7-I-1.2 in units of 4.5 kpc. If the cloud is roughly spherical and has uh ~ lO^cm"^, 
then objects will have an additional ~ 5 magnitudes of extinction with respect to objects on 
the near side. Thus, a natural explanation for the geometry of G70.7+1.2 is that the bright 
Be-star is moving into the near side of the cloud creating a reflection nebula, while on the 
far side, a heavily extincted B-star /pulsar binary is moving into the cloud creating a bow 
shock and powering the nebula (Figure [3]). 

One remaining puzzle is the origin and impact of the hot gas powering the diffuse X- 
ray emission. Figure [T] shows that the radio/optical and diffuse X-ray morphologies are 
substantially different, and the peak of the diffuse X-ray emission is separated ~ 20"from 
the center of the radio/optical emission (which contains the Be-star and X-ray binary). 
Thus, it is apparent that this hot gas does not play an important dynamic nor, given the its 
luminosity is ~ 10^^ erg s~^, energetic role in shaping G70.7-I-1.2. Two viable explanations 
for the origin of the hot gas are, given the quasi-shell like morphology, that it is the result 
of a previous explosive event that the X-ray binary is overtaking or it may be unassociated 
with the G70.7-I-1.2. In any case, the origin of this plasma — either related or unrelated to 
G70.7-I-1.2 — is unknown. 

The definitive proof of the proposed model (Figure [3]) would be the discovery of a 
pulsar associated with G70.7-I-1.2. A search for pulsations with the Green Bank Telescope 
at 2.2 GHz is underway. If a pulsar is found, G70.7-I-1.2 will be an important laboratory for 
studying plasma processes taking place between the pulsar/B-star wind and the interaction 
of that mixture with the cold molecular gas. 
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Table 1. X-ray Spectral Fits. 



Model 


Nh 


V/kBT 


Flux 






(10^2 cm-2) 


(keV) 


(10~^^ erg cm~^ s~^) 




(1) 


(2) 


(3) 


(4) 


(5) 


Diffuse Emission 


Power-law 


0.701°:}^ 




311?? 


18.2/13 


Rayiiioud- Siiiit li Plasma 


i.oi;-? 


0.7112:- 


331^ 


9.75/13 


CXO J200423.4+333907 


Power-law 


-L-'J-l.O 


1.81- 


1 7+5.6 
-L- ' -1.7 




Power-law (P = 2.0) 




(2.0) 


J^-"-0.4 




Power-law (A^^ = 1 x 10^^) 


(1.0) 


1 q+0.5 
-'-•'J-0.4 







Note. — All errors are 68% confidence levels. Values in parentheses are held fixed. (1) - 
Absorbed spectral model. (2) - Best-fit column density. (3) - Measured photon index for 
power-law models and ksT for the Raymond-Smith plasma. (4) The unabsorbed flux in the 
0.5-2.5 keV band for the diffuse emission and 2-10 keV band for CXO J200423.4+333907. 
(5) - The value of for diffuse emission models and the number of degrees of freedom, 
V. This column is not apphcable to CXO J200423.4-I-333907 since the spectral fitting was 
performed with Cash statistics. 




Fig. 1. — H-band image of G70.7+1.2 with contours of the adaptively smoothed X-ray 
emission (black hnes) and the Chandra 99% confidence (0'.'49) error circle (blue circle). The 
X-ray contours are logarithmically spaced between 10% and 90% of the peak emission. 
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Fig. 2. — Color magnitude diagram using data from iBessell fc BrettI (1l988l ). Filled circles 
show the IR counterpart for Ay = 5.5 and Ay = 10.5 at a distance of 4.5 kpc . 
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Fig. 3. — Diagram of the geometry of G70.7+1.2. The newly discovered X-ray binary moves 
into the far side of the molecular material powering the radio/NIR/optical nebula, whereas 
the Be-star creates a reflection nebula on the near-side. See the text for details. 



